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(57) Abstract 

Imaging tumors using diffuse light. An imaging system includes a source (35) of diffuse light for generating oscillatory diffuse 
photon density waves to illuminate an object (10). a detector (60) for detecting diffuse photon density waves interacting with the object, 
and a computer (130) interfaced with the detector for processing data corresponding to the photon density waves detected to determine at 
least a position of the object. In one embodiment, the turbid medium (20) and the object have associated therewith at least one diffusion 
coefficient and the diffuse photon density waves which illuminate the object refract around the object as a result of their interaction with it, 
thereby producing a distorted wavefrom that allows the computer to construct an image of the object. In another embodiment, a fluorescent 
object (100) produces re-radiated diffuse photon density waves which allow the object to be imaged. 
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OBJECT IMAGING USING DIFFUSE LIGHT 

Field of the Invention 

This invention relates generally to imaging of 
objects. More specifically, this invention relates to methods 
5 and apparatus for imaging objects using diffuse light. 

Background of the Invention 

Techniques for imaging objects have been used for 
nearly a century in the medical arts for diagnosing and 
understanding the myriad diseases and maladies that afflict 

10 the human body. Imaging techniques have also found use in 

such diverse fields as radio astronomy, sonar, radar and other 
fields which require information about an object which is not 
readily visible to the naked eye and therefore not easily 
examined. Medical imaging techniques include, for example, X- 

15 ray imaging, positron emission tomography (PET) , ultrasound 
imaging and the well known magnetic resonance imaging (MRI) 
In all of the imaging techniques mentioned above, 
narrow band frequency radiation illuminates the object of 
interest to produce reflected or emitted radiation which is 

20 then gathered from the object by a detector. The reflected or 
emitted radiation is then processed by an imaging algorithm to 
obtain useful information about the object. 

In medical applications, the use of ionizing 
radiation in imaging, for example with X-rays, involves 
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significant health risks to a patient when the patient is 
exposed to the radiation for prolonged periods of time or in 
multiple imaging schemes. Furthermore, certain of these 
imaging techniques undesirably involve the use of invasive 
5 procedures which are both costly and painful. Yet other 
techniques such as MRI do not yield consistently useful 
clinical results. 

There has thus arisen in the medical imaging art an 
interest in developing non-invasive/ safe and relatively fast 

10 techniques which can take advantage of the natural scattering 
of visible and infrared light through media containing objects 
to be imaged. Techniques using diffuse light could be used 
in conjunction with other imaging schemes such as X-ray 
imaging or MRI to produce highly useful clinical images for 

15 diagnostic purposes. 

Much of the progress in imaging with diffusive light 
has focused on ballistic techniques using lasers. With these 
techniques, an intense pulsed laser illuminates a sample. By 
time gating photons that have been scattered only a few times 

20 and rejecting all other photons, the optical absorption of the 
medium and objects found therein can be mapped. This 
technique works best when the allowed time window is short and 
the photons deviate the least from their "ballistic" 
trajectory. Unfortunately, the transmittal intensity of 

2 5 unscattered photons diminishes exponentially with increasing 

sample thickness. 

Because of the limitations of ballistic imaging, it 
is difficult to obtain high quality images of relatively thick 
objects with low power lasers. Examples of ballistic imaging 
30 techniques are disclosed in K. M. Yoo, F. Lie and R. R. 

Alfano, Optics Letters, Vol. 16, p. 1068 (1991), and in D. A. 
Benaron and D. K. Stevenson, Science, Vol. 259, p. 1463 
(1993) . 

A second technique for imaging using diffuse light 

3 5 is optical phase modulation. Phase modulation techniques have 

permitted the location of single absorbers using low power, 
continuous wavelength lasers . In accordance with these 
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techniques, an amplitude modulated source creates photon 
density waves that acquire anomalous phase shifts due to the 
absorber. For the case of a single absorber, the distortions 
are readily interpreted; however for a more complicated object 
5 a general analysis is required. 

An example of imaging with diffuse light is 
disclosed in U.S. Patent No. 5,119,815, Chance where scattered 
light was applied to a biological imaging application. The 
Chance patent discloses a technique for solving the diffusion 

10 equation for a homogeneous medium to obtain the overall 

optical absorption characteristics. This was possible for the 
homogeneous medium because the long time limit of the 
logarithmic derivative of the detected intensity yields the 
absorption characteristics directly. Thus the absorption 

15 characteristics for uniform structures may be obtained with 
the methods and apparatus disclosed in the Chance patent. 

Still other attempts to image with diffuse light are 
disclosed in U.S. Patent No. 5,070,455, Singer et al . In the 
Singer et al . system, light intensities are measured at many 

20 sensor positions (pixels) , initial values of absorption or 

scattering coefficients are assigned at each pixel, and then a 
new set of intensities at each pixel is calculated. The 
calculated intensities are compared to the real intensities, 
and the intensity differences are used to generate a 

25 subsequent interaction of absorption or scattering values for 
each pixel . 

The methods described in Singer et al . usually 
require many iterations since the absorption or scattering 
values may not converge rapidly. Furthermore, the Singer et 
30 al. system utilizes cumbersome Monte-Carlo statistical 

techniques which consume large amounts of processing time 
without guaranteeing computational success. Singer et al . ' s 
methods may also produce false local minima providing 
misleading results for the absorption characteristics. 
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Thus prior imaging techniques using diffuse light 
for scattering fail to solve a long-felt need in the art for 
robust imaging techniques which can produce reliable images in 
biological systems. Solution of the aforementioned problems 
5 has heretofore eluded the medical imaging art. 

Summary of the Invention 

The aforementioned deficiencies in the imaging art 
are overcome by methods and apparatus provided in accordance 
with the present invention which provide imaging of objects in 

10 turbid media using diffuse light. 

In a preferred embodiment of the invention an 
imaging system comprises source means for generating 
oscillatory diffuse photon density waves to illuminate the 
object, detection means for detecting diffuse photon density 

15 waves produced as a result of the diffuse photon density waves 
interacting with the object, and processing means interfaced 
with the detection means for processing data corresponding to 
the photon density waves detected by the detection means to 
determine at least a position of the object. The turbid 

2 0 medium and the object have associated therewith at least one 

diffusion coefficient and the diffuse photon density waves 
which illuminate the object diffract around or refract through 
the object as a result of their interaction with it, thereby 
producing a distorted wavefront such that after the detection 
25 means detects the distorted wavefront the processing means 

determines the diffusion coefficient of the turbid medium and 
the object. More preferably, the processing means constructs 
phase- and amplitude contours corresponding to propagation of 
the distorted wavefront and further determines at least the 

3 0 position of the object from the phase and amplitude contours, 

thereby imaging the object. 

In yet a further preferred embodiment of imaging 
systems provided in accordance with the invention, display 
means are interfaced with the processing means for displaying 
3 5 the image of the object, the source means comprises at least 
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one laser, and the detection means comprises an optical fiber 
interfaced with a photomultiplier tube. 

Further aspects of the invention provide imaging of 
a . fluorescent object such that the diffuse photon density 
5 waves having a first wavelength cause the object to fluoresce, 
thereby producing re-radiated diffuse photon density waves 
having a second wavelength such that after the detection means 
detects the re-radiated diffuse photon density waves, the 
processing means can image the object, 

10 In the embodiment of the invention where 

fluorescent, re-radiated diffuse photon density waves are 
detected, the source means preferably comprises a plurality of 
lasers oriented around the object which alternately irradiate 
the object with the diffuse photon density waves of the first 

15 wavelength to cause the object to fluoresce. The detection 
means comprises an optical fiber that is placed in proximity 
to the object and a photomultiplier tube interfaced to the 
optical fiber. The imaging system further comprises switch 
means interfaced with each of the plurality of lasers for 

2 0 alternately and sequentially turning on and off each laser, 
and radio frequency driving means interfaced through the 
switch means with the lasers for driving the lasers to produce 
the diffuse photon density waves of the first wavelength. 

Alternatively, the imaging system comprises a 

2 5 plurality of lasers each having a spatial location with 

respect to the object. Each laser is more preferably 
modulated at all times during imaging at a different frequency 
in a frequency range around a specified frequency, thereby 
producing a power spectrum associated with each spatial 

3 0 location around the object. Analysis means are provided 

interfaced with the detection means and the processing means 
for analyzing the power spectrums associated with each spatial 
location to determine the position of the object. 

In still a further preferred embodiment of imaging 
3 5 systems which take advantage of fluorescent, re-radiated 
diffuse photon density waves, the source means comprises a 
phased-array . The phased-array , preferably comprises at least 
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two lasers that are substantially one hundred and eighty 
degrees out of phase with each other, thereby producing the 
diffuse photon density waves having the first wavelength which 
interfere destructively to produce an amplitude null line and 
5 a substantially one hundred and eighty degree phase shift 

across the null line equidistant from the lasers. By scanning 
the null line with the phased-array the processing means can 
produce an image of the object in the turbid medium. 

Systems and methods provided in accordance with the 

10 present invention will provide efficient imaging of tumors and 
other maladies that effect human tissue. These systems will 
also prove to be much more economical to build when compared 
to prior imaging systems, since they will not require the 
complex machinery that is associated with prior imaging 

15 systems such as MR I . Furthermore, since the methods and 

apparatus provided in accordance with the present invention 
utilize diffuse light for imaging, more reliable images of 
tumors will be presented to the medical diagnostician or 
clinician so that cancerous tumors and other inhomogeneities 

20 in the tissue will be detected earlier and more readily. This 
holds the promise of saving lives and reducing the overall 
costs of medical care. 

The features, objects and advantages of the 
invention will be understood by those with skill in the art 

2 5 from the following detailed description of the preferred 

embodiments thereof when read in conjunction with the drawings 
'which are first described briefly below. 

Brief Description of the Drawings 

Figure 1 is a schematic representation of an imaging 

3 0 system using diffuse light wherein scattering experiments with 

diffuse photon density waves occurs. 

Figure 2 depicts measured constant contours of 
diffuse photon density waves propagating through a turbid 
medium. 

35 Figure 3 depicts measured constant phase contours of 

diffuse photon density waves propagating through turbid media 
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and refracting across a plane boundary between two turbid 
media with different light diffusion coefficients. 

Figure 4 is a graph of phase contours of diffuse 
photon density waves propagating through turbid media after 
5 refracting across a cylindrical boundary separating two turbid 
media wherein lensing effects are observed. 

Figure 5 is a schematic representation of a "time- 
sharing" system for imaging diffuse photon density waves re- 
radiated from a fluorescent object. 
10 Figures 6A is a graph of phase contours of the 

disturbance produced by a fiber point source in a turbid 
medium. ...... 

Figure 6B is a graph of the absorption and emission 
characteristics of a fluorescent dye contained in a spherical 
15 shell within the turbid medium. 

Figure 6C is a graph of measured constant amplitude 
phase contours for diffuse photon density waves emitted from a 
source into a turbid medium, and re-radiated diffuse photon 
density waves. 

20 Figure 7 is a schematic representation of a 

"frequency-encoded" system for imaging diffuse photon density 

waves re-radiated from a fluorescent object. 

Figures 8A and 8B are schematic representations of 

"phased-array " scanning systems for imaging re-radiated 
25 diffuse photon density waves and graphs of the object's 

position in the turbid medium and a null line associated with 

the phased-array. 

Detailed Description of Preferred Embodiments 

Propagating disturbances are produced in a dense, 
3 0 turbid medium containing objects when amplitude modulated 
light sources are introduced into the turbid medium. In 
biological and medical imaging applications, human tissue, 
such as breast tissue containing tumors, is such a turbid 
medium containing objects to be imaged. In accordance with 
35 the present invention, an oscillatory light source introduces 
diffuse photon density waves (hereinafter referred to as 
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"DPDW") in the turbid medium. While the results described 
herein have been achieved with respect to several experimental 
apparatus to be described hereafter, those with skill in the 
art will immediately recognize that the techniques and 
5 apparatus disclosed are readily applicable to imaging human 
body tissue which is infected with tumors or other maladies 
that will appear as inhomogeneities in the tissue which is a 

turbid medium. 

The inventors of the subject matter herein claimed 

10 and disclosed have discovered that DPDW can be used in at 
least two ways to perform imaging. In a first preferred 
embodiment, DPDW are introduced into a turbid medium and 
refracted through and diffracted around objects in the medium, 
thereby producing a distorted wavefront which can then be 

15 analyzed to yield useful information concerning at least a 
position and size of the object around which the DPDW have 
been refracted. 

In a further preferred embodiment, the inventors 
have determined that by introducing fluorescence to the 

20 object, photons in the re-radiated DPDW will be emitted from 
the object at a shifted wavelength with respect to the 
original DPDW and the re-radiated DPDW then analyzed to 
determine at least the position of the fluorescent objects in 
the turbid medium. This is particularly useful for imaging of 

25 tumors since in prior, non-diffuse light imaging techniques 
fluorescent dyes which respond to various forms of radiation 
have typically been introduced into tumors, thereby yielding 
information about the location and size and nature of the 
tumor under examination. 

3 0 Thus, imaging tumors and objects in turbid 

media with diffuse light will in the future produce 
significant and effective images for diagnostic and clinical 
purposes . 

I. THE BASIC THEORY OF DPDW PROPAGATION IN A 
3 5 TURBID MEDIUM 

The basic mechanics of wave propagation in a 

turbid medium have been explored by the inventors when 



NSDCX3D: <WO 9512132A1> 



i 



WO 95/12132 * 




PCT/US94/12486 



- 9 - 

performing imaging in accordance with the present invention. 
DPDW are scalar, over-damped, traveling waves of light energy 
density, denoted U(r,t) . They propagate through any medium in 
which the transport of light energy density, U, is governed by 
5 the "diffusion equation, " which is 




= DV 2 U 



where D is the "diffusion coefficient" for the medium. This 
10 diffusion equation holds f or ■ a non-absorbing medium. -Some 
examples of optically turbid media include dense suspensions 
of micrometer-sized spheres, human tissue, paints, foams, and 
Intralipid (a mixture of water, soybean oil and egg yolk) . 
The introduction of amplitude modulated light into a turbid 
15 medium produces a macroscopic ripple of brightness that is 
microscopically composed of individual photons undergoing 
random walks. The disturbances arise whenever the diffusive 
system is driven by an oscillating source. ;> 

The oscillatory part of the solution for an in- 
20 finite, non-absorbent, homogeneous turbid medium in the 

presence of a point source located at the origin follows the 
form: 

U(r,t) = (A/Dr) (exp{ -kr } ) (exp{ i (kr - cox}) 

where A is a constant, r is the radial distance from the 
25 origin, D is the diffusion coefficient for light in the 
medium, CO is the source modulation frequency, and k = 
(co/2D) 1/2 . Although the wave is very rapidly attenuated, it 
has a well-defined wavelength, amplitude and phase at all 
points. Interestingly, the wavelength can be altered by 
30 modifying D or co. 

When absorption is present, a similar solution 
to Equation (1) can be obtained, but the real and imaginary 
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parts of the wave vector k are different and depend explicitly 
on the sample absorption length (as well as the photon random 
walk step for inverse scattering factor) . The macroscopic 
disturbance obeys a Helmholtz equation, and therefore has many 
5 properties that are normally associated with conventional 
electromagnetic radiation. Thus, DPDW display refraction, 
diffraction, and interference properties similar to those 
observed with conventional electromagnetic radiation 
propagating through media. By examining the energy contours 
10 produced as a result of the interaction of the DPDW with the 
objects in the medium according to the techniques of the 
present invention, the objects can' be imaged to obtain useful 
clinical data. 

II. IMAGING USING REFRACTED DPDW 

15 In a preferred embodiment, an apparatus for 

imaging objects in a turbid medium is shown schematically in 
Figure 1. The object 10 may contain smaller objects such as 
tumors which must be characterized for clinical purposes. For 
experimental purposes, the dense turbid medium imaged is 

2 0 Intralipid, a polydisperse suspension of particles having an 
average diameter of -0.4 Jim, but a relatively wide range of 
sizes (i.e., from -0.1 to -1.1 Jim). By changing the solution 
concentration of the Intralipid, it is possible to vary the 
light diffusion coefficient, D, of the medium. The photon 

2 5 transport mean free path/* was about 0.2 cm in a 0.5% 

concentrated solution. 

To perform imaging of the smaller objects in 
the Intralipid, a large fish tank 20 (30 cm x 30 cm x 60 cm) 
is filled with this material. When the absorption is very 
30 small, the suspensions in the Intralipid are dilute, and 

therefore the diffusion coefficient is inversely proportional 
to the Intralipid concentration. 

In a preferred embodiment, a source fiber 3 0 
and detector fiber 40 (both -4 mm in diameter) are placed in 

3 5 proximity to object 10. The source light 3 0 is preferably 
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derived from a 3-mW diode laser 35 operating at about 816 run. 
The diode laser 3 5 is amplitude modulated at 2 00 MHz by driver 
50, and the position of the source fiber 3 0 is fixed. The 
detector fiber 40 is in detecting proximity to the object 10 
5 and is further connected to a photomultiplier tube 60 on its 
other end. In order to facilitate phase and amplitude 
measurements, both the reference from the driver 50 and the 
detected signal are down- converted to 25 kHz by heterodyning 
with a second oscillator 70 to 200.025 MHz. Mixers 80 

10 preferably combine the frequencies to produce the 2 5 kHz 

signals shown. The low-frequency signals are then measured 
• using a lock-in amplifier 90. 

In a more preferred embodiment, the- phase shift 
(and ac amplitude) of the detected light is measured with 

15 respect to the source 30 at each point on a 0.5-cm square 
planar grid throughout the sample. Constant phase contours 
are then determined by linear interpolation of the grid. data. 
The sensitivity of the apparatus shown in Figure 1 is. about 
10 5 . Since the signal amplitude decays by >e' 2n in one 

20 wavelength, the range of the apparatus is limited to slightly 
more than one wavelength. Nevertheless, it is possible -to 
clearly distinguish the essential physical phenomena of the 
DPDW scattering through the object 10 with this apparatus. 

The results obtained by imaging with the system 

25 of Figure 1 demonstrate that a "diffusional index of 

refraction" exists for DPDW, and that it is possible to 
manipulate the diffusional index of refraction by controlling 
the photon diffusion coefficients (D) of adjacent turbid 
media. This is of considerable importance in biological 

30 systems, where the natural curvature of organs such as the 
brain, heart, or kidney, together with changes of scattering 
and absorption as in the grey-white matter transition of the 
brain, can lead to significant modifications of the 
trajectories of diffuse photons. 

35 Results for a -0.5% concentrated homogeneous 

turbid medium are illustrated in Figure 2. Constant phase 
contours are shown at 20° intervals about the source. The 
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contours are approximately circular, and their radii can 
therefore be extrapolated back to the source 30. In the 
insert of Figure 2, the phase shift and the quantity 
lnl rC7 ac (rt)l as a function of radial distance from the source 
5 are plotted. From these measurements it is possible to 
determine that the wavelength of the DPDW is 11.2 cm, the 
photon transport mean free path is - 0.2 cm, and the photon 
absorption length is -52 cm in -0.5% Intralipid at 22°C. In 
this case, the photon absorption can be attributed almost 

10 entirely to water. 

Referring to Figures 3 and 4, refraction of the 
DPDW is shown in three ways: (1) when a plane, acrylic 
partition or (2) a cylindrical acrylic partition is placed in 
the tank to simulate a boundary, and (3) when no partition is 

15 in place so that the medium is homogeneous. In Figure 3, 
constant phase contours exist at about every 20°. When the 
plane boundary 45 is introduced, the lower medium 55 has a 
concentration, c lt - 1.0% and light diffusion coefficient D x . 
The upper medium 65 has a concentration, c u , « 0.25% and light 

20 diffusion coefficient D u . The contours below boundary 45, 
shown at the 4 cm position on the y-axis, are just the 
homogeneous media contours without refraction. These contours 
are obtained before the partition is introduced. The contours 
above the boundary 45 are derived from the DPDW that are 

2 5 transmitted into the less concentrated medium. 

Theory predicts that the wavelength in the less 
dense medium, 7^, is 14.8 cm, and should be greater than the 
wavelength of the DPDW in the incident medium, X lt which is 
8.17 cm. The ratio of the two wavelengths should equal the 

3 0 ratio of the diffusional indices of refraction of the two 

media. As is predicted, X u = X 2 (D x /D u )~ in - X x (c 2 /c u ) 1/2 . 

Theory also predicts that the apparent source 
position S if as viewed from within the upper medium 65, should 
be shifted from the real source position, S 0 = 4.0^+0.2 cm, by 
35 a factor = 0.55. This is verifiable from the data of 

Figure 3 since, using the radii from the full contour plots, 
the apparent source position is shifted from 4.0 +. 0.2 to 2.0 
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± 0.25 cm . 

In accordance with the invention, Figure 3 
explicitly demonstrates Snell's law for DPDW. This can be 
seen by following the ray from S a to the point A at the 
5 boundary 45, and then into the upper medium 65. The ray in 
the lower medium 55 makes an angle 0^ = 14° with respect to 
the surface normal. The upper ray is constructed in the 
standard way between the apparent source position S if through 
the point A on boundary 45, and into the medium 65 above the 

10 boundary. The upper ray is perpendicular to the circular wave 
fronts in the less dense medium 55, and makes an angle 9 C = 
26.6° with respect to the boundary normal. Therefore, it can 
be seen graphically that sin^/sinG,. = 0.54 ~ X 1 /X ut so that 
Snell's law accurately describes the propagation of DPDW 

15 across boundary 45. 

The inventors have also determined that by 
using a circular boundary (shown generally at 75 in Figure 4) 
to separate the two turbid media, the curvature of the DPDW 
can be altered in analogy with a simple lens in optics*. 

20 Referring to Figure 4, two semi-infinite media are separated 
by curved boundary 75, and the medium 55 on the right is more 
concentrated. The constant phase contours of the transmitted 
wave exhibit a shorter wavelength, and are clearly converging 
toward some image point to the right of the boundary. The 

25 medium 65 on the left (Xj) has an Intralipid concentration of 
-0.1%, and the medium 55 on the right {X T ) has a concentration 
of -1.6%. 

The wavelength ratio is measured to be X r /^ 2 = 3.8+. 
0.3. The curved surface 75 has a radius S = 9.0 + 0.4 cm. 

30 The position of source 30 is S 0 = 9.4 ± 0.3 cm. The image 
position is determined to be S x = 12 + 2 cm. This result 
deviates somewhat from the well-known paraxial result from 
geometrical optics for imaging by a spherical refracting 
surface. The deviation is primarily a result of spherical 

3 5 aberration. However, in accordance with the invention, the 
curvature of the wave fronts is reversed after traversing the 
circular boundary 75. 



1NSOOCID:<WO 95121 32A1> 



WO 95/12132 




PCT/US94/12486 



- 14 - 

The results shown in Figures 2, 3 and 4 which 
are obtained from the apparatus of Figure 1 show that it is 
possible to exert substantial control over the transport of 
diffuse light in dense random media. The inventors have 
5 clearly demonstrated that the index of refraction of DPDW in 
such a medium depends on the photon diffusion coefficient or 
random walk step of the photons in the medium. This allows 
imaging of inhomogeneities in the medium by examining the 
DPDW, a result that has not heretofore been achieved in the 
10 art. 

It should be noted that when deriving Equation 
(1), the differential from of Fick's law and photon flux 
conservation principles were used. In a preferred embodiment 
of imaging by examining the scattering and diffraction 

15 refraction of DPDW, the effects of absorption of the DPDW can 
be ignored, and it can be assumed that the time it takes for 
light to travel a single random walk step is much shorter than 
the modulation period. In this case, the oscillatory part of 
the light energy density ^(r) obeys the Helmholtz equation, 

20 i.e., (V 2 + k 2 ) C7 w (r) =0 . The only significant difference of 

DPDW propagation in comparison to conventional wave phenomena 
is that k 1 = i(co/D), and therefore k is complex. 

The spatial part of L7 ac (r,t) in Equation (1) is 
simply the Green's function solution of the Helmholtz equation 

25 with the appropriate k. Therefore, some of the basic theorems 
that apply to solutions of the Helmholtz equation will apply 
to DPDW propagation in a turbid "medium. For example, a 
Kirchoff integral can be constructed for these waves using the 
Green's function solution. This provides a formal method by 

3 0 which to calculate the wave amplitude and phase at various 

distances from a "diffracting" aperture as has been discussed 
in this Part II. To the extent that the Kirchoff integral 
embodies the basic Huygens-Fresnel principle, contributions of 
different elements of a scattering surface arising from 

3 5 damped, spherical point sources will be observed. This also 
implies that the focusing of DPDW will have the same 
limitations due to diffraction as in the case of light 
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propagating in a standard optical medium. Thus, imaging of 
tumor-like inhomogenei ties in tissue can be accomplished in 
accordance with the present invention by examining the 
refractive, diffractive and scattering properties of DPDW 
5 incident to the tissue. 

III. IMAGING BY EXAMINING RE-RADIATED DPDW 

The inventors of the subject matter herein 
claimed and disclosed have also discovered that examining re- 
radiated DPDW from a fluorescent inhomogeneity in a turbid 
10 medium provides methods of imaging inhomogeneities in the 
* medium. ' Many types of tumors which occur in the human body 
comprise inhomogeneities which can be made to fluoresce after 
being irradiated with DPDW, and therefore, the inventors have 
determined that re-radiated DPDW provide an excellent means of 
15 locating tumors. 

In order to observe re-radiated DPDW and image 
an object, the model biological material Intralipid is used. 
An amplitude modulated 200 MHz or 50 MHz laser diode ~3mW, 780 
nm is fiber-coupled into the medium/ and another optical fiber 
20 is used to detect diffuse photons as a function of position 
within the medium. Using the standard heterodyne techniques 
discussed earlier, the phase and amplitude of the DPDW in the 
medium can be observed. 

A preferred apparatus for imaging inhomo- 
2 5 geneities in turbid media by detecting re-radiated DPDW is 

shown schematically in Figure 5. The absorber /radiator shown 
at 100 preferably comprises a shell filled with Intralipid at 
a concentration of 0.4mg/L Indocyanine green dye which is less 
than one tenth of the concentration commonly used in human 
30 subjects to test hepatic function. In this embodiment, 

multiple sources 110 and a single detector 120 are used to 
determine the center of the fluorescent object 100. 

Computer processor 13 0 preferably outputs a 
control signal to switch 140 which controls the sequential 
35 activation of each of sources 110. A detector 150 which most 
preferably comprises a photomultiplier tube is interfaced both 
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to detector fiber 120 and computer processor 130. In a 
further preferred embodiment, detector 150 is a Hamamatsu R928 
or Rl645u photomultiplier tube which also comprises a high 
voltage power supply that ensures adequate gain to computer 
5 processor 13 0 so the computer processor 13 0 can process the 
data received from detector fiber 120. A display device 160, 
preferably a CRT screen, is interfaced to computer processor 
130 and outputs an image 170 of object 100. 

To obtain image 170, a fluorescent object 100 

10 is irradiated with multiple sources 110 of DPDW. It is then 
preferred to measure the amplitude (or phase) of the re- 
radiated light. The partial amplitude resulting solely from 
source i, is dependent on the i th source-object separation, 
the quantum efficiency of the dye in the object 100, and the 

15 object 100-detector 120 separation. Therefore: 

exp(-kl r i - ~r Q I ) e xp(-kTr Q - r d I ) 

I Uj ~ £, — 

20 I r ± - ~r 0 I I r 0 - r d I 

(2) 



where £, is the quantum efficiency of the dye, r t is the 

25 position of the i th source, r 0 is the position of the object 
center, r d is the detector position, and k (k' ) is the wave 
vector magnitude of the DPDW at 780 run (83 0 nm) . 

3 0 The individual sources 110 are separately 

turned on and off, and the re-radiated amplitude for each 
object separation is measured. Since the source positions r ± , 

35 and the detector position r d are known, it is possible to 

estimate the object's position by finding the value of r G that 
gives the best agreement with the measured ratio I /I Ujl . 
Generally, three sources 110 are necessary to localize the 

40 object in two dimensions. However, it is more preferable to 
use four or more sources to improve the signal-to-noise ratio 
of the imaging system. Using four sources 110 as shown in 
Figure 5, it is possible to localize the center of the 1 cm 
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spherical object 100 to within 0.4 cm. This two-dimensional 
localization is easily extended to three dimensions, as will 
be recognized by those with skill in the art. 

Referring to Figure 6A, constant phase contours 
5 of the disturbance produced by the fiber point source located 
at the origin of the system are illustrated. In a more 
preferred embodiment, re-radiated DPDW are obtained by filling 
a spherical glass shell with the absorbing dye Indocyanine 
green, and then illuminating the sphere with DPDW in the 
10 Intralipid solution. The dye is preferably chosen to absorb 
radiation at the source wavelength of 7 80 nm, and very soon 
thereafter re-radiate photons at a red-shifted energy, 83 0 nm. 
Because the dye has a lifetime of less than 1 nsec compared to 
the 5nsec period of the source, the re-radiated energy is in 
15 the form of a DPDW at the red-shifted energy. 

The absorption and emission characteristics of 
the dye are shown in Figure 6B wherein the inset is the 
chemical formula of the dye material. The Intralipid solution 
surrounding the obstacle has a concentration of 0.1% giving a 
20 source diffuse photon density wavelength of -18 cm. A point 
source at the origin generates the incident DPDW. Constant 
amplitude contours of the incident wave in the presence of the 
obstacle are shown in Figure 6C (dashed lines) . In a more 
.preferred embodiment, two spectral filters centered on 83 0 nm 
25 (Schott glass filters, RG830) are provided to enable the 
incident and re-radiated DPDW to be separated. 

In Figure 6C, the measured constant amplitude 
contours of the wave at 830 nm are shown as solid lines, and 
the measured incident amplitude contours at 7 80 nm are shown 
30 as dashed lines. The dashed lines demonstrate the DPDW 
character of the re-radiated waves. The re-radiated wave 
originates from within the absorbing object. As can be seen 
from the contours, the DPDW wavelength at 83 0 nm is somewhat 
longer than the DPDW wavelength at 780 nm. This is a function 
35 of the relatively larger diffusion coefficient for 830 nm 

light in Intralipid. Thus, a type of fluorescence of DPDW has 
occurred and the inhomogeneity is converted into a source of 
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secondary DPDW. 



An alternative embodiment to the " time- sharing 



apparatus" for imaging using re-radiated DPDW described by 
Figure 5 is a "frequency encoded apparatus" shown in Figure 7. 
5 In this embodiment, each source 12 0 is modulated by a slightly 
different modulation frequency, fi, (i.e., f A = 2 00.00 + k 
(.01) MHz) and is kept on at all times. In this way, a 
modulation power spectrum is associated with each spatial 
location in the sample. By measuring the power spectrum of 

10 the re-radiated light as a function of modulation frequency 
with a spectrum analyzer 180, the position of object 100 can 
be determined using essentially the same analysis as described 
above with respect to the apparatus of Figure 5. This type of 
imaging is similar to MRI in that a frequency spectrum emitted 

15 by an object is converted to the spatial location of the 
object . 



fluorescent objects that re-radiate DPDW is shown in Figures 
8A and 8B. The apparatus of Figures 8A and 8B is 

2 0 qualitatively different than the time-sharing apparatus and 
frequency encoded apparatus of Figures 5 and 7 respectively, 
since it uses a scanning phased-array 190 and a single 
detector 200. Preferably, the phased-array 190 comprises two 
sources 85 and 95 that are substantially 180° out of phase 

2 5 with each other and that emit DPDW which interfere 

destructively to produce an amplitude null and a sharp 180° 
phase shift across a curve that describes this family of 
points called the "null line," shown schematically at 210. 



3 0 and then moving an absorbing object 220 from one side of the 
null line 210 to the other, the object 220 will preferentially 
absorb light from the nearest source, and therefore distort 
the null line. When the absorber 220 is also a re-radiator, 
the complimentary effect is seen, that is, the object re- 

35 radiates more light derived from the closest source. In both 
measurements, the phase of the detected DPDW will undergo a 
180° shift as the absorber crosses the original, undisturbed 



A third, alternative embodiment for imaging 



By placing a detector 2 00 on the null line 210, 
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null line 210. These effects are demonstrated by the graphs 
of Figure 8A. The re-radiated light displays a sharper phase 
transition, and a deeper amplitude null than that of the 
incident light, and the re-radiated wave exhibits a 
5 complimentary change as discussed above. 

Interestingly, in the embodiment of Figure 8A, 
the phase shift of the re-radiated light is always the same, 
independent of detector position. Thus, it is possible to 
hold the absorber 220 and the detector 200 stationary, and 

10 scan the null line 210 which will produce essentially the same 
results. Scanning the null line 210 can be achieved by 
mechanically translating the- two sources 85 and 95 together. 

. Referring to Figure 8B this effect is 
demonstrated for re-radiated light. A sharp phase shift at 

15 the location of the re-radiator 220 and null line 210 is 

detected. The phase transition of the incident light occurs 
near the position of detector 200. With knowledge of the 
position of the null line 210 as a function of time, a one- 
dimensional localization of the re-radiator can be achieved. 

20 By performing three scans down three perpendicular axes,, 
three-dimensional localization of re-radiator 220 can be : 
achieved. It is envisioned that the phased-array apparatus of 
Figures 8A and 8B will be simple to build and implement. 
Furthermore, imaging apparatus described herein will generally 

25 be much more economical to build than other imaging systems 
currently in use, such as MRI systems, PET systems, and CAT 
* scanners . 

Therefore, it has been demonstrated that imaging of 
inhomogeneities in turbid media using DPDW is both possible 

3 0 and practical with the methods and apparatus disclosed and 

claimed herein. By studying the scattering and distortion of 
the DPDW by inhomogeneities found therein, the inhomogeneities 
can be located and their sizes determined for clinical 
purposes. Alternatively, a fluorescent inhomogeneity can be 

3 5 examined by analyzing the re-radiated DPDW emitted therefrom. 
In accordance with the invention, it is also possible to 
combine the dual techniques of examining the refractive, 
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diffractive, scattering, and fluorescent properties of the 
inhomogeneities to image the objects. 

The techniques and apparatus of the present 
invention provide particularly powerful clinical tools for 
5 imaging, locating and sizing tumors in human tissue. 

Furthermore, the methods and apparatus claimed and disclosed 
herein should prove to be both economical and efficient 
imaging tools for clinical diagnosis. 

There have thus been described certain preferred 

10 embodiments of methods and apparatus for imaging objects using 
diffuse light provided in accordance with the present 
invention. While preferred embodiments have been described 
and disclosed, it will be recognized by those with skill in 
the art that modifications are within the true spirit and 

15 scope of the invention. The appended claims are intended to 
cover all such modifications. 
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CLAIMS 

What is claimed is: 

1. A system for imaging an object using diffuse 

light comprising: 
5 source means for generating oscillatory diffuse 

photon density waves to illuminate the object; 

detection means for detecting diffuse photon 
density waves produced as a result of the diffuse 
photon density waves interacting with the object; 
10 and 

processing means interfaced' with the detection 
means for processing data corresponding to the 
photon density waves detected by the detection means 
to determine at least a position of the object. 

15 2 . The system recited in claim 1 wherein the 

turbid medium and the object have associated therewith at 
least one diffusion coefficient and the diffuse photon density 
waves which illuminate the object diffracts around or refracts 
through the object as a result of their interaction with it, 

2 0 thereby producing a distorted wavefront such that after the 

detection means detects the distorted wavefront the processing 
means determines the diffusion coefficient of the turbid 
medium and the object. 

3 . The system recited in claim 2 wherein the 
25 processing means constructs phase and amplitude contours 

corresponding to propagation of the distorted wavefront and 
further determines at least the position of the object from 
the phase and amplitude contours, thereby imaging the object. 

4. The system recited in claim 3 further 

3 0 comprising display means interfaced with the processing means 

for displaying the position of the object. 

5. The system recited in claim 4 wherein the 
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source means comprises at least one laser and the detection 
means comprises an optical fiber interfaced with a 
photomultiplier tube . 

6. The system recited in claim 1 wherein the 
5 object is fluorescent and the diffuse photon density waves 
which have a first wavelength cause the object to fluoresce, 
thereby producing re-radiated diffuse photon density waves 
having a second wavelength such that after the detection means 
detects the re-radiated diffuse photon density waves, the 
10 processing means can image the object. 



7. The system recited in claim 6 wherein the 
source means comprises a plurality of lasers oriented around 
the object which alternately irradiate the object with the 
diffuse photon density waves of the first wavelength to cause 

15 the object to fluoresce. 

8. The system recited in claim 7 wherein the 
detection means comprises an optical fiber that is placed in 
proximity to the object and a photomultiplier tube interfaced 
to the optical fiber. 

9. The system recited in claim 8 further 
comprising switch means interfaced with each of the plurality 
of lasers for alternately and sequentially turning on and off 
each laser and radio frequency driving means interfaced 
through the switch means with the lasers for driving the 
lasers to produce the diffuse photon density waves of the 
first wavelength. 

10. The system recited in claim 9 further 
comprising display means interfaced with the processing means 
for displaying the position of the object produced by the 

3 0 processing means. 



MSDOCID: <WO 9512132A1> 



WO 95/12132 



PCTAJS94/12486 



- 23 - 

11. The system recited in claim 6 wherein the 
source means comprises a plurality of lasers each having a 
spatial location with respect to the object and each laser 
being modulated at all times during imaging at a different 

5 frequency in a frequency range around a specified frequency, 
thereby producing a power spectrum associated with each 
spatial location around the object. 

12. The system recited in claim 11 further 
comprising analysis means interfaced with the detection means 

10 and the processing means for analyzing the power spectrums 
associated with each spatial location to determine the 
position of the. object. 

13 . The system recited in claim 12 further 
comprising display means interfaced with the processing means 

15 for displaying the image of the object. 

14. The system recited in claim 6 wherein the 
source means comprises a phased- array . 

15. The system recited in claim 14 wherein the 
phased-array comprises at least two lasers that are 

2 0 substantially one hundred and eighty degrees out of phase with 
each other, thereby producing the diffuse photon density waves 
having the first wavelength which interfere destructively to 
produce an amplitude null line and a one hundred and eighty 
degree phase shift across the null line. 

25 16. The system recited in claim 15 wherein the 

phased-array scans the null line so that the processing means 
can produce the image. 
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17. A method of imaging an object in a turbid 



medium using diffuse light comprising the steps of: 



illuminating the object with oscillatory 
diffuse photon density waves; 



5 



detecting diffuse density waves which are 
produced as a result of the diffuse photon density 
waves interacting with the object; and 



10 



determining at least the position of the object 
in the turbid medium by analyzing the detected 
diffuse photon density waves. 



18. The method recited in claim 17 wherein the 



turbid medium and the object have associated therewith at 
least one diffusion coefficient and the diffuse photon density 
waves which illuminate the object refract around the object as 
15 a result of their interaction with it, thereby producing a 
distorted wavefront which is analyzed to determine the 
diffusion coefficient of the turbid medium and the object. 



determining step comprises the steps of constructing phase 
20 contours corresponding to propagation of the distorted 

wavefront and determining at least the position of the object 
from the phase contours, thereby imaging the object. 

20. The method recited in claim 19 further 
comprising the step of displaying the image of the object. 

25 21. The method recited in claim 17 wherein the 

object is fluorescent and the diffuse photon density waves 
which have a first wavelength cause the object to fluoresce, 
thereby producing re-radiated diffuse photon density waves 
having a second wavelength such that detecting the re-radiated 

30 diffuse photon density waves the object to be imaged. 



19 . 



The method recited in claim 18 wherein the 
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22. The method recited in claim 21 wherein the 
illuminating step comprises the step of alternately 
irradiating the object with at least two sources of diffuse 
photon density waves of the first wavelength to cause the 

5 object to fluoresce. 

23. The method recited in claim 21 wherein the 
illuminating step comprises the step of continuously 
illuminating the object with at least two sources of diffuse 
photon density waves each having a spatial location with 

10 respect to the object, wherein each source has associated 

therewith a' frequency in a frequency range around a specified 
frequency, thereby producing a power spectrum associated with 
each spatial location around the object. 

24. The method recited in claim 21 wherein the 
15 illuminating step comprises the steps of irradiating the 

object with a phased-array and scanning an amplitude null line 
to produce the image of the object. 

25. A system for imaging an object in a turbid 
medium comprising: 

2 0 source means for illuminating the object with 

oscillatory diffuse photon density waves of a first 
specified wavelength, whereby the object will 
fluoresce re-radiate diffuse photon density waves of 
a second wavelength after being illuminated with the 
25 oscillatory diffuse photon density waves of the 

first specified wavelength; 

detection means for detecting the re-radiate 
diffuse photon density waves of the second 
wavelength; and 

3 0 processing means interfaced with the detection 

means for processing data corresponding to the re- 
radiated diffuse photon density waves of the second 
wavelength to determine at least the position of the 
object in the turbid medium. 
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26. The system recited in claim 2 5 wherein the 
source means comprises a plurality of lasers oriented around 
the object which alternately irradiated the object with the 
diffuse photon density waves of the first wavelength to cause 

5 the object to fluoresce. 

27. The system recited in claim 26 wherein the 
detection means comprises an optical fiber that is placed in 
proximity to the object and a photomultiplier tube interfaced 
to the optical fiber. 

10 28. The system recited in claim 27 further 

comprising switch means interfaced with each of the plurality 
of lasers for alternately and sequentially turning on and off 
each laser and radio frequency driving means interfaced 
through the switch means with the lasers for driving the 

15 lasers to produce the diffuse photon density waves of the 
first wavelength. 



29. The system recited in claim 28 further 
comprising display means interfaced with the processing means 
for displaying the image of the object produced by the 
2 0 processing means. 



30. The system recited in claim 25 wherein the 
source means comprises a plurality of lasers each having a 
spatial location with respect to the object and each laser is 
modulated at all times during imaging at a different frequency 

25 in a frequency range around a specified frequency, thereby 
producing a power spectrum associated with each spatial 
location around the object . 

31. The system recited in claim 30 further 
comprising analysis means interfaced with the detection means 

3 0 and the processing means for analyzing the power spectrums 
associated with each spatial location to determine the 
position of the object. 
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32. The system recited in claim 31 further 
comprising display means interfaced with the processing means 
for displaying the image of the object. 

33. The system recited in claim 25 wherein the 
5 source means comprises a phased-array. 



phased-array comprises at least two lasers that are 
substantially one hundred and eighty degrees out of phase with 
each other, thereby producing the diffuse photon density waves 
10 having the first wavelength ' which interfere destructively to 
produce an amplitude null line and a one hundred and eighty 
degree phase shift at points equidistant from the lasers. 

35. The system recited in claim 34 wherein the 
phased-array scans the null line so that the processing means 

15 can produce the image. 

36. A method of imaging an object in a turbid 
medium comprising the steps of: 

illuminating the object with oscillatory 
diffuse photon density waves of a first specified 
2 0 wavelength; 



34 . 



The system recited in claim 33 wherein the 



allowing the object to fluoresce, thereby re- 
radiating a diffuse photon density waves having a 
second wavelength; 



25 



detecting the re-radiated diffuse photon 
density waves having the second wavelength; and 



analyzing the photon density waves having the 
second wavelength to determine at least the position 
of the object in the turbid medium. 
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37. The method recited in claim 36 wherein the 
illuminating step comprises the step of alternately 
irradiating the object with at least two sources of diffuse 
photon density waves of the first wavelength to cause the 

5 object to fluoresce. 

38. The method recited in claim 3 6 wherein the 
illuminating step comprises the step of continuously 
illuminating the object with at least two sources of diffuse 
photon density waves each having a spatial location with 

10 respect to the object, wherein each source has associated 

therewith a frequency in a frequency range around a specified 
frequency, thereby producing a power spectrum associated with 
each spatial location around the object. 

39. The method recited in claim 36 wherein the 
15 illuminating step comprises the steps of irradiating the 

object with a phased-array and scanning an amplitude null line 
to produce the image of the object. 

40. A system for imaging an object in a turbid 
medium wherein the turbid medium and object have a diffusion 

20 coefficient associated therewith comprising: 

source means for illuminating the object with 
oscillatory diffuse photon density waves, whereby a 
distorted wavefront of diffuse photon density waves 
is created by refraction of the oscillatory diffuse 
25 photon density waves refracting around the object; 

detection means for detecting the distorted 
wavefront of diffuse photon density waves refracted 
around the object; and 

processing means interfaced with the detection 
3 0 means for processing data corresponding to the 

distorted wavefront of diffuse photon density waves 
to determine at least the position of the object in 
the turbid medium. 
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41. The system recited in claim 40 wherein the 
processing means constructs phase contours corresponding to 
propagation of the distorted wavefront and further determines 
at least the position of the object from the phase contours, 

5 thereby imaging the object. 

42. The system recited in claim 41 further 
comprising display means interfaced with the processing means 
for displaying the image of the object. 

43 . The system recited in claim 42 wherein the 
10 source means comprises at least one laser and the detection 
means comprises an optical fiber interfaced with a 
photomultiplier tube . 

44. A method of imaging an object in a turbid 
medium wherein the turbid medium and the object have a 

15 diffusion coefficient associated therewith comprising the 
steps of: 

illuminating the object with oscillatory 
diffuse photon density waves; 

allowing the diffuse photon density waves" to 

2 0 refract around the object throughout the turbid 

medium thereby producing a distorted wavefront of 
diffuse photon density waves; 

detecting the distorted wavefront of diffuse 
photon density waves; and 
25 analyzing the distorted wavefront of diffuse 

photon density waves to determine at least the 
position of the object in the turbid medium. 

45. The method recited in claim 44 wherein the 
determining step comprises the steps of constructing phase 

3 0 contours corresponding to propagation of the distorted 

wavefront and determining at least the position of the object 
from the phase contours, thereby imaging the object. 
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46. The method recited in claim 45 further 
comprising the step of displaying the image of the object. 
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